Decay of Persistent Spin Helix due to the Spin Relaxation at Boundaries 
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We study electron spin relaxation in one-dimensional structures of finite length in the presence of 
Bychkov-Rashba spin-orbit coupling and boundary spin relaxation. Using a spin kinetic equation 
approach, we formulate boundary conditions for the case of a partial spin polarization loss at the 
boundaries. These boundary conditions are used to derive corresponding boundary conditions for 
spin drift-diffusion equation. The later is solved analytically for the case of relaxation of a homoge- 
neous spin polarization in ID finite length structures. It is found that the spin relaxation consists 
of three stages (in some cases, two) - an initial D'yakonov-Perel' relaxation is followed by spin 
helix formation and its subsequent decay. Analytical expressions for the decay time are found. We 
support our analytical results by results of Monte Carlo simulations. 
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Dynamics of electron spin polarization in semiconduc- 
tor structures has attracted a lot of attention recently 
in the context of spintronicJ^^, which is playing a fun- 
damental role in the novel technological developments 
based on different effects in this scientific area. More- 
over, the ability to understand and predict the dynamics 
of electron spins in semiconductors is also important for 
the area of two terminal electronic devices with mem ory, 
so-called memristive device^®. In some of thenPSl, the 
electron spin degree of freedom defines their internal state 
and, consequently, is responsible for their time-dependent 
memory response. 

It has been shown by us recentlj'sHIII that the system 
boundaries significantly modify the dynamics of process. 
For example, we have demonstrated^^ that in finite size 
2D systems, the spin polarization density decays much 
slower than in the bulk and the exponential spin po- 
larization decay rate is defined by both the system size 
and strength of spin-orbit interaction. In finite length 
wire^ and channels-^ (oriented in a specific direction), 
changes in the electron spin relaxation are even more pro- 
nounced: instead of relaxing to zero, the homogeneous 
electron spin polarization relaxes into a persistent spin 
polarization structure known as the spin helij^^^ - a 
spin polarization configuration in which the direction of 
spin polarization density rotates along the wire. 

In real experimental situations, the spin helix config- 
uration can not exist infinitely long. Here, we assume 
that the main decay mechanism is due to spin relaxation 
at system boundaries. Indeed, local strong random elec- 
tric fields in the vicinity of boundaries result in a ran- 
dom spin-orbit interaction influencing the electron spin 
degree of freedom. It is thus important to develop a the- 
ory and model spin relaxation in constrained geometries 
taking into account the boundary spin relaxation and un- 
derstand how the boundary spin relaxation changes the 
overall character of electron spin relaxation in the entire 
system. 

In this paper, we use both spin kinetic!^ and diffu- 
gjQjjiSHial equations to investigate the dynamics of elec- 



tron spin polarization in semiconductor wires of finite 
length. Specifically, we consider dynamics of spin relax- 
ation in one-dimensional (ID) finite length systems with 
Bychkov-Rashbc^ spin-orbit interaction and boundary 
spin relaxation. Since it is easier to incorporate the 
boundary spin relaxation into the boundary conditions 
for spin kinetic equations^^, below, we use the spin ki- 
netic equation approach first. Next, based on bound- 
ary conditions for the spin kinetic equations, we derive 
boundary conditions for spin diffusion equation which 
is easier to solve. Finally, we obtain an exact solution 
for the problem of spin relaxation in finite length wire. 
Our analytical studies are complemented by semiclassi- 
cal Monte Carlo simulations of spin dynamics^-'^ giving 
an additional insight into the problem. 

A kinetic description of electron spin polarization in 
one-dimensional wires can be giverfl^ in terms of vectors 
S"'' and S^, which are the spin polarizations of electrons 
moving along the wire in the positive (with momentum 
p = mvex), and negative (p = —mvex) x-directions with 
the average velocity v = 1/t, where I is the mean free 
path and r is the momentum relaxation time. The ki- 
netic equation can be written as a system of two vector 
equation^ 

which take into account electron spin precessions (first 
term in the r.h.s. of Eqs. ([T]), ^) induced by 
Bychkov-RashbcPl spin-orbit interaction and bulk scat- 
tering events (second term in the r.h.s. of Eqs. ([I]), ([2|). 
Here, = 2ap/h is the angular velocity (directed along 
y-axis) of spin precessions, a is the spin-orbit coupling 
constant, p ~ mv is the average momentum, ey is the 
unit vector in y-direction. Additionally, it is convenient 
to introduce the parameter r] = Q,/v = 2amh~^, which 
gives the spin precession angle per unit length. Note that 
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the direction of spin precessions in Eqs. Q, ([2]) is differ- 
ent for left- and right-moving electrons. 

Eqs. ([!]), ([2]) incorporate the well-known D'yakonov- 
Perel' spin relaxation mechanisn i^^^^ - l Indeed, trajecto- 
ries of spin-polarized electrons are randomized bulk scat- 
tering events described by second term in the r.h.s. of 
Eqs. ([1]), ([2]). Correspondingly, the direction of spin 
rotation becomes fluctuating what causes average spin 
relaxation (dephasing). The incomplete spin relaxation 
in finite-length wirea^J can be explained by the existence 
of a maximum electron spin precession angle defined by 
the system size. 

In this work, we study the dynamics of spin polar- 
ization in a wire of the length L, — L/2 < a; < i/2, 
in the presence of spin relaxation at the boundaries 
r ~ [x — ±L/2]. The boundary conditions, which take 
into account the relaxation of spin polarization in elas- 
tic scatterings of electrons from the boundaries, can be 
formulated as 



[S- = 7S + ] 
[S+ = U=_L/2, 



(3) 
(4) 



where 7,0<7<1, isa phenomenological dimensionless 
papameter characterising the boundary relaxation rate. 
If 7 = 1, then the spin polarization of scattered electrons 
is conserved (no boundary spin relaxation). The case 7 — 
corresponds to the total relaxation of spin polarization 
at the boundary, whereas for < 7 < 1 we deal with a 
situation of a partial boundary relaxation. 

Combining Eqs. (Tj), ([2| we obtain a single equation 
for the total spin polarization S = S+ + S~ 



1 dS 

T~dt 



dx'^ 



9S 

— 2ilvey X — 
dx 



(S - ^^e^) = 0, (5) 



where Sy is the y-component of S. What are the bound- 
ary conditions for S? These can be derived reformulating 
the boundary conditions given by Eqs. (3), Q in terms 
of the function S only. Subtracting Eq. ( 2 1 from Eq. ([!]) 
we find 



^(S+ 
dt ^ 



as 

dx 



f2e„ X S + - 

T 



0.(6) 



Now, let us consider, for example, the boundary x = L/2. 
From Eq. ^ we see that &t x = L/2, S+ = S/(l + 7) 
and S~ — 78/(1 -|- 7). Substituting these boundary val- 
ues into Eq. (6]) we obtain the corresponding bound- 
ary condition that can be generally formulated for both 
boundaries as 



dx 1 + 7 \ ot T 



0.(7) 



Introducing a complex polarization 
S — Sx + iSz , 



(8) 



it is straightforward to rewrite Eq. ([5]) and boundary 
conditions ([7| in a more compact form 
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S^O, (9) 
-0. (10) 



x=±L/2 



The diffusion limit of Eqs. ^ and ( 10 1 is realized when 
L ^ I, rjl 1, and t ^ t. In this case we can neglect 
d'^S/dt^ compared to T'^dS/dt in Eq. ^ and dS/dt 
compared to S/ r in Eq. ( 10 1 , because all relevant spin 



relaxation times are much longer than the momentum 
relaxation time r. As a result, we obtain 



dS^ 
~dt 
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x=±L/2 



The initial condition for Eq. ( 11 ) is 



S{x,t^ 0) = So{x). 



(11) 
(12) 

(13) 



We can exclude rotations of spin polarization vector 
(defined by Eq. (Isl)) that are still present in Eq. (11) 



and Eq. (12 I by introducing a complex field u via 



u{x, t) 



'S{x,t). 



(14) 



Note that the rotation transformation ( |14[ ) is a particular 
case of more general transformations that can be use d to 
remove spin precessiorP^. It can be showrP that Eq. (11) 



and Eq. (12 1 transform into the ordinary heat equation 



du 
dt 

with boundary condition 



D 



d'^u 
dx"^ ' 



7 u 



du ^\ 
dx 1 + 7 Z 



(15) 



(16) 



x=±L/2 



where D — P jr. Note that in the diffusion approxima- 
tion j/-component of spin polarization, Sy, satisfies the 
same Eq. (flSj) and Eq. ([T6|). 



We find the general solution of Eq. ( 15 ) with the 



boundary conditions ( 16 ) using the standard method of 



x=±L/2 



separation of variables. The straightforward application 
of this method leads to the following expression for the 
complex spin polarization 



It is important to keep in mind that Sy is not coupled 
to any other component of spin polarization (see Eqs. 
([5]), ([7|). Consequently, we can safely take out Sy from 
our consideration selecting Sy{x,t = 0) = 0. 



S = 



plTjX 



(17) 
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FIG. 1: Relaxation of homogeneous spin polarization ini- 
tially polarized in ^-direction, (a) and (b) show z- and x- 
components of spin polarization density at several moments 
of time. The wire length L = lOOZ, r]l = 0.1545, 7 = 0.97. 
These parameters result in k = 1.31, and ^1 = 0.78. 



where /i„ and ^,1 are positive roots of 

tan/i„ = — K/.t„, and cot^„ = k^„. 



(18) 



andK = 2(l+7);/[(l-7)L]. 

Thinking about relaxation of an initial spin polariza- 
tion profile, it is evident that an 'overall' spin relax- 
ation rate is determined by the smallest of the roots of 
Eqs. (18 1, which is the root ^1 satisfying the inequal- 
ity < ^1 < 7r/2 (note that 7r/2 < ^1 < tt). We 
can find explicit expressions for the spin relaxation time 
Tr = L^/{4:D£^f) in the limiting cases of small and large 
K. Specifically, when k <C 1 (the limit of strong boundary 
spin relaxation), 



TT 

= and Tr 



(1 + 2k). (19) 



In the opposite limit, when k 3> 1 (the limit of weak 
boundary spin relaxation), we obtain 

6k -1 J (1+7)(1 + (3k)-1)Lt 



OK 2 



2(1 -7)Z 



(20) 
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FIG. 2: Temporal evolution of the spin polarization density 
Sz at the wire's center for the same parameter values as in 
Fig. [1] The fitting curves are obtained employing numerical 
values of tdp and Tr calculated for these parameter values. 



There are two additional time scale in the problem - 
the D'yakonov-Perel' spin relaxation time, t^p = t j (Irj)'^ 
and the time of persistent spin helix formation in the 
absence of boundary spin scattering, Th- As these time 
scales are hidden in the set of relaxation times describing 
dynamics of the general solution (Eq. (17)), we refer 
to our prior publicatiorP where the expression for t/j is 
given, Th = L"^ /{tt^D). It is interesting that th can not 
exceed t^, namely, > t/j. Consequently, when Tr ^ 
T/i, the spin dynamics can be considered as a three-stage 
process: an initial D'yakonov-Perel' spin relaxation (i) is 
followed by spin helix formation (ii), which is followed 
by its decay {in). Otherwise, when and are close, 
the decay stage significantly overlaps with the spin helix 
formation and thus the stages {ii) and {Hi) can not be 
well separated. 

In the case of three-stage dynamics, the spin helix am- 
plitude (at the end of the second stage) is defined by 
the relation between t^p and t/j times. In particular, if 
"Th 3> Tp,p (or, equivalently, Lrj 3> tt), then the spin he- 
lix amplitude is much smaller than the amplitude of the 
initial spin polarization. In the opposite regime, when 
Th Tdp (or Lrj < tt), the spin helix amplitude is com- 
parable to the initial one. 

As an application of the above theory, let us consider 
a specific problem - the problem of relaxation of initially 
homogeneous spin polarization pointing in z direction. In 
this case, the initial condition simply reads S{x,0) = iSa 
(see Eq. ([s])). The coefficients a„ and 6„ are obtained 
by substitution this initial condition into Eq. ( |17[ ), its 
multiplication by an appropriate sine or cosine function 
and subsequent integration over x. This results in 
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FIG. 3: Comparison of z-component of spin polarization 
calculated analytically (straight lines) and numerically using 
Monte Carlo simulations (dots) for L = lOOZ, r]l = 0.01545, 
7 = 0.8. 



.2^„ sin^„(cos(0.577i) — O.STyLrt sin(0.5r?L)) 
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(22) 

Fig. [l] demonstrates dynamics of initially homogeneous 
spin polarization plotted using Eqs. (17 1, (21), (22 1. We 



emphasize that the initial relaxation at the center is of 
D'yakonov-Perel' type as, e.g., Sz and Sx in the central 
region (see t = lOOr curves) are flat. The spin helix 
configuration is clearly seen dX t ^ IOOOt. The subse- 
quent slower evolution of spin polarization is caused by 
the boundary spin scattering. The boundary spin scat- 
tering results in smaller amplitudes of spin polarization 
oscillations closer to the ends and a larger amplitude at 
the wire's center. 

Three stages of spin polarization dynamics, however, 
are better visualized plotting the magnitude of Sz at 
the wire's center as a function of time. Fig. [2] shows 
Sz{Q,t) for two selected values of boundary spin scatter- 
ing coefficient 7. Clearly, both curves contain two in- 
tervals of exponential relaxation that are seen on this 
logarithmic plot as straight lines. The initial interval of 
D'yakonov-Perel' spin relaxation (0 < t < 200t) is not 
influenced by the boundary spin relaxation. The later, 
however, determines the character of spin relaxation at 
long times 750r < t. The interval of spin helix formation 



(200t < i < 750t)) is located between two regions of 
exponential evolution mentioned above. 

In order to obtain an additional insight on spin relax- 
ation with boundary spin scattering, we have performed 
Monte Carlo simulations using an approach described in 
Refs. |25] and EH This approach is based on a semiclas- 
sical description of electron space motion and quantum- 
mechanical description of spin dynamics. All main parts 
of the Monte Carlo simulation algorithm are described in 
Refs. US] and [m and will not be repeated here. The only 
novel feature of the code is the boundary spin scatter- 
ing mechanism implemented in the following way. When 
an electron scatters from a boundary, the direction of its 
spin is inverted with a probability p = (1 — 7)/2. It is 
not difficult to notice that this part of the algorithm is 
equivalent to the boundary conditions (Eqs. (|3|, ([4|) for 
the spin kinetic equation. 

There is an excellent agreement between our analytical 
results and the numerical Monte Carlo simulations. For 
example, let us consider dynamics of spin relaxation in 
a finite length wire with boundary spin scattering when 
the spin precession angle per mean free path is relatively 
small. We also assume that at the initial moment of 
time i = 0, the spin polarization is homogeneous and 
points in z direction. For this situation. Fig. [3] presents 
a comparison of Sz at different moments of time plotted 
using Eqs. (17), (21), (22) (smooth curves,) and Monte 
Carlo simulations (dotsJT^ We see that there is an ex- 
cellent agreement between the analytical and numerical 
results. 

In conclusion, we have developed a theory of spin relax- 
ation in wires accounting for boundary spin relaxation. 
For both spin kinetic and diffusion equations appropriate 
boundary conditions have been derived. Based on this 
theory, we predict the existence of three (in some cases, 
however, two) stages of spin dynamics consisting of an 
initial D'yakonov-Perel' relaxation followed by spin he- 
lix formation and its subsequent decay. Experimentally, 
parameters of boundary spin relaxation can be extracted 
from both the long-time spin helix decay rate and spin 
helix shape distortion. 
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